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ABSTRACT 

Since a majority of young low-mass stars are members of multiple systems, 
the study of their stellar and disk configurations is crucial to our understand- 
ing of both star and planet formation processes. Here we present near-infrared 
adaptive optics observations of the young multiple star system VWCha. The 
previously known 0'.'7 binary is clearly resolved already in our raw J and K band 
images. We report the discovery of a new, faint companion to the secondary, 
at an apparent separation of only O'.'l or 16 AU. Our high-resolution photometric 
observations also make it possible to measure the J — K colors of each of the three 
components individually. We detect an infrared excess in the primary, consistent 
with theoretical models of a circumprimary disk. Analytical and numerical cal- 
culations of orbital stability show that VWCha may be a stable triple system. 
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Using models for the age and total mass of the secondary pair, we estimate the 
orbital period to be 74 years. Thus, follow-up astrometric observations might 
yield direct dynamical masses within a few years, and constrain evolutionary 
models of low-mass stars. Our results demonstrate that adaptive optics imaging 
in conjunction with deconvolution techniques is a powerful tool for probing close 
multiple systems. 

Subject headings: binaries: close — circumstellar matter — infrared: stars - 
stars: formation — stars: low-mass, brown dwarfs — stars: pre-main sequence 

1. Introduction 

Most stars are members of binary or multiple systems (Batten 1973; Halbwachs 1983). 
Yet the formation of multiple stellar systems is still poorly understood. Currently, the most 
promising mechanism appears to be fragmentation during the collapse of dense molecular 
cloud cores (Bodenheimer et al. 2000). Recent surveys of pre-main-sequence (PMS) stars 
in certain star-forming regions have measured a binary frequency about twice that in the 
solar neighborhood (see Mathieu et al. (2000) for a review and caveats). Whether this is an 
evolutionary or environmental effect is open to debate. In any case, higher-order multiple 
systems appear to be relatively rare in both young and old stellar populations (see, e.g., 
Ghez et al. (1997) and references therein). 

High-resolution observations of PMS multiple systems help determine the stellar and 
disk configurations of young objects, leading to better constraints on models of multiple star 
formation and early evolution. We have undertaken an adaptive optics imaging program 
of known PMS close binary systems to search for additional companions and to detect 
circumstellar or circumbinary disks. In this Letter, we present our observations of one such 
system, VWCha in the ChamaeleonI star-forming cloud at a distance of ~160pc (Whittet et 
al. 1997). The VW Cha system has previously been reported to consist of several components; 
Schwartz (1991) listed Sz23 located 16" away as a possible wide companion; Brandner et al. 
(1996) discovered a closer companion at 0"7 separation; Ghez et al. (1997) observed these 
components and reported an additional companion at 2"7 separation from the primary. 
Brandner & Zinnecker (1997, BZ97) obtained spatially resolved spectra of the close binary. 

Here, we present the discovery of a new close companion to the secondary of VW Cha 
and an infrared excess from the primary in the system, probably due to a dusty disk. 
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2. Observations 

We observed the VWCha system on 11-13 March 2000 with the European Southern 
Observatory 3.6-meter telescope at La Silla, Chile, using the adaptive optics near infrared 
system (ADONIS) and the SHARPII+ camera. The SHARPII+ camera is based on a 
256x256 NICMOS III array. We used a plate scale of 0.035 arcsec/pixel, Nyquist-sampling 
the diffraction limited point spread function (PSF) and giving an effective field of view 
of 8''5x8'.'5. J band (1.25 /xm) observations of VWCha were conducted during the first 
night and K band (2.18 /im) during the second. At J we obtained 60 frames, each with an 
integration time of 10 s, chopping the source out of the field of view for a 10 s integration on 
sky between each frame. The science target itself was used as the wavefront sensor of the 
adaptive optics (AO) system. The seeing during the VWCha observations the first night, 
as reported by the La Silla seeing monitor, fluctuated in the V band (0.55 /im) between 0'.'8 
and 1"4. The K band observations of the second night were obtained in a similar manner 
with 10 frames of 30 s integration time each, with the difference that the sky frames were 
obtained after the series of on-source frames and not in between. The seeing in V this time 
was stable at about 1"0. During both nights most frames showed a diffraction limited core 
while the Strehl ratios fluctuated between 0.1 and 0.3. 

After each series of observations of VW Cha, an identical series of observations with the 
same AO correction parameters were obtained of a PSF calibrator star, for post processing 
purposes. We used the F5V star SAO 256804 with V magnitude 8.50. 

3. Reductions and Analysis 

The basic data reductions were done in a standard way by subtracting sky frames from 
source frames and then dividing by a flat field obtained on sky during dusk. The resulting 
frames were then not stacked together but instead processed with the myopic deconvolution 
algorithm ID AC 4 . 

From the PSF calibrator star observations, we obtained an initial guess of the PSF 
during the science data acquisition. Since the PSF calibrator star data were not obtained 
simultaneously with the science data, the initial guess PSF is then modified by IDAC it- 
eratively to better match the different frames of science data. A series of PSF estimates 
corresponding to the science object series is thus produced along with the deconvolved im- 



4 Information about IDAC and the package itself may be found on the web pages of ESO, 
http:/ /www. eso.org/ 
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age. The deconvolution procedure using IDAC took about two weeks of processing time on 
an UltraSparc workstation. We refer to Christou et al. (1999) for details on IDAC. 

We performed astrometry of the VW Cha system using the deconvolved images in both 
J and K. The error in separation was estimated by checking the consistency of measured 
separations of other multiple star systems observed several times during the two nights. 
The error in position angle was then estimated from the error in separation, by assuming 
positional errors to be isotropic. 

For flux calibration, we used observations of 10 other SAO stars distributed over the sky 
and obtained during the two nights in both J and K. Their known B and V magnitudes and 
spectral classes were used to derive their J and K magnitudes using standard colors (Bessell 
& Brett 1988) and an estimate of the color excess due to extinction using Table 1 of Mathis 
(1990), where we assumed an optical extinction ratio Ry = 3.1. We then fitted a relation 
to the observed count rates and derived magnitudes, estimating the absolute photometric 
errors from the statistical error of the fit. By checking the consistency of the measured flux 
of stars observed several times during the two nights we derived the relative photometric 
errors. 



4. Results and Discussion 

4.1. Multiplicity of VW Cha 

In our observations, we are able to detect binaries at separations between 0'.'08 and 4", 
where the lower limit corresponds to the diffraction limit of the 3.6-m telescope in J, and 
the upper limit is due to the field of view. The inner 0'.'7 binary is already resolved 

in the raw data in both J and K. After deconvolution, the secondary is split up into a tight 
O'.'IO binary in itself; see Fig. 1. Our results are summarized in Table 1. 

The reported 2"7 companion was not seen in any of our images. Based on the background 
noise levels, our 5a detection limits are J = 20 and K = 17, well below K « 9.4 ±0.6 derived 
from Table 2 in Ghez et al. (1997). The cause for a drop in brightness by a factor of 10 3 is 
hard to imagine and we suggest that either the star has drifted out of our field of view, or 
that the detection was spurious. 
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4.2. Implausibility of chance alignment 

Comparing our observations with measurements from 1992 (Brandner et al. 1996), we 
find no significant change (less than 0"05) of the relative positions between the A and B+C 
components. Since VWCha has a proper motion of about / .'22yr _1 (Frink et al. 1998), the 
system as a whole has moved at least 2" over the past 10 years. Thus the proper motions 
of the components are highly correlated, corresponding to a relative velocity of less than 
5kms _1 at the distance of VWCha. We therefore conclude that the VWCha system is not 
due to chance alignment but a physical triple star. 

4.3. Stellar parameters 

4-3.1. Corrections for extinction 

From observations in the visual (BZ97) and in the infrared (Table 1) it follows that 
the primary increasingly dominates the emission at increasingly longer wavelengths. In- 
terpolating between the flux ratios in different wavebands, we aportion the observed total 
/-magnitude = 10.60 (Gauvin & Strom 1992) as 11.03 to A and 11.73 to B+C. The resulting 
observed colors, (/ — J) a = 1-82 and (/ — J)b+c — 1-79, can then be compared to the intrin- 
sic colors of K5/K7 and K7 stars (BZ97) for logg of 3.5 and solar metallicity of the model 
atmospheres by Allard et al. (2000, AHS00), viz. (I - J) = 0.80 and 0.92 respectively. 
The color excess of the relatively 'blue star' B+C is probably mainly due to foreground dust 
extinction, hence Aj — Aj = 0.87 which yields Ay =3.03 mag, where we have used the ex- 
tinction curve for the Chal cloud (G. Olofsson, private communication: e.g. Ay = 2.90 Aj, 
Aj c = 1.83 Aj, Ak = 0.37 Aj). The assumption that this value of Ay also applies to the 
primary would imply excess emission from component A, at the level of about 0.1 mag at J 
and 0.9 mag at K. An active circumprimary disk would be a natural source of such excess 
radiation. 



4-3.2. Stellar luminosities and masses 

For the extinction and excess emission corrected J-magnitude of the primary, J 0j a — 
8.16, the absolute bolometric magnitude can be obtained from [AHS00] as Mb i,A = 3.6. For 
B and C, this value is 5.0 and 5.3, respectively. Consequently, the bolometric luminosities are 
L A = 2.9 L , L B = 0.79 L Q and L c = 0.59 L Q . The 'blue' components B and C contribute 
to the total flux mainly at shorter wavelengths. The neglect of correcting for extinction is 
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expected, therefore, to underestimate the total luminosity of the system. This is indeed the 
case, as the calorimetric luminosity of L ca \ = 3L by Prusti et al. (1992) constitutes merely 
70% of Lboi = 4.3 L . A natural explanation of this deficit would be that only a fraction of 
the stellar photons are reprocessed to the infrared, while the rest escapes non-isotropically 
due to the presence of a circumstellar disk. 

The empirically determined parameters of the individual stellar components of VW Cha 
are consistent with those of the theoretical evolution models (for single stars) 5 by Palla & 
Stahler (1999) for masses of A, B and C of 1.0 M , 0.40 M Q and 0.35 M , respectively, at the 
common age of 0.4 x 10 6 yr (Table 1). This assumes a systematic shift of AlogT e ff = —0.04, 
not unreasonable. 

BZ97 found the Ha-line to be in emission toward both A and B+C. If dominated by 
disk accretion processes, some infrared excess could be expected also from the fainter pair 
and should thus be corrected for. However, the Ha flux is larger from A (BZ97) and the 
mass accretion rate, M, is presumably also larger for the primary. Estimates of M, which 
are in reasonable agreement with the K-excess, are in the range of a few times 10" 8 M yr _1 
(Gahm et al. 1995) to 10 _7 M Q yr _1 (Hartmann et al. 1998), but certainly much less than 
the large rates (10^ 6 to 10" 4 M o yr- 1 ) obtained by Johns-Krull et al. (2000) in the far-UV, 
a spectral region which is highly susceptible to extinction corrections. 

VWCha was detected by ISO at 6.7 /xm and 14.3 /xm (Persi et al. 2000), but not by IRAS 
at 60 /xm and 100 /xm (Prusti et al. 1992), suggesting that only small amounts of dust exist 
in the system which are at temperatures significantly below 50 K. For the stellar properties 
and mass accretion rate determined above, an outer circumprimary disk radius < 70 AU 
could then be inferred from standard theory of steady state, thin-disk accretion (Frank et al. 
1985), using a temperature profile of the form T^^r) oc r~ 0A . This minimum temperature 
would also be consistent with the scenario of a tidally truncated disk around A of size of 
at most twice the BC binary separation, advocated in the next section (the separation of 
0'.'2 corresponds to 30 AU, where T disk ~ 70 K). Further, the observed excess in J would be 
consistent with the presence of hot dust and thus of only a small central hole. However, 
as illustrated by the interpretation of the observed weakness of the 10 /xm silicate feature in 
VW Cha (Natta et al. 2000) some caution seems to be called for, when applying simplified 
disk theories developed for single stars to close multiples. 

It is interesting to note that the secondary in the TTau binary system recently has 
been resolved into a binary in itself (Koresko 2000). The projected relative distances for 
A:(B+C) and B:C are about 100 AU and 10 AU respectively, making the configuration very 



5 The validity of these models for binaries has recently been assessed by Palla & Stahler (2001). 
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similar to VWCha. In addition, evidence for a circumprimary disk around TTau with a 
radius ~40 AU has been presented by Akeson et al. (1998). The major difference is the much 
higher luminosity of the T Tau system, ~30 L (Cohen & Kuhi 1979). An apparent difference 
is also the extremely red appearence of TTauS. This has previously been interpreted as the 
(unresolved) secondary being in a considerably earlier stage of evolution than the primary. 
However, later results suggest that the IR-excess may be due to foreground extinction (see 
discussion in Akeson et al. (1998)). 



4.4. Triple star dynamics: Stability estimates 

We consider the possibility that VWCha is a physical triple, bound and stable for 
~1 Myr (upper limit on age). Components B and C form a tight pair B:C at the projected 
separation of 16 AU, while A is found at a 7 times larger projected separation of 106 AU. 
We first show analytically that a nearly planar 3-body system similar to VW Cha, with low 
eccentricities, is forever stable against the disintegration of the close binary by the third 
massive body. 

We utilize energetic constraints on the motion of a 3-body system, similar to the familiar 
zero-velocity curves in the restricted 3-body problem where the Jacobi integral entails the 
so-called Hill stability, involving no exchange of components and no 'ionization'. A small 
binary B:C will never be pulled apart by component A, if (Marchal and Bozis (1982), Marchal 
(1988)) 

(-EL 2 )>(IU 2 ) L1 , (1) 

where E is the total mechanical energy of the triple system, L its total angular momentum, 
/ the semi-moment of inertia (one-half times the sum of the squares of radius vectors times 
masses), and U the total potential energy of the triple. Subscript LI denotes the value at the 
collinear "Lagrange point" between the two most massive bodies. Stability does not depend 
on either the total mass or the distance units. E.g., tilting a given system with respect to 
the plane of the sky merely rescales the inferred periods but not the shapes of orbits. If 
projection effects are not extreme, the orbital periods are close to those in a planar system 
with semi-major axes of a\ = 110 AU and a 2 = 16 AU for the A:(B+C) and B:C pairs, i.e. 
on the order of 872 yr and 74 yr. 

For the masses given in Table 1, we found that right hand side of the stability criterion 
equals (IU 2 )li = 0.0152GM t 2 ot (where M tot = 1.75 M is the total mass). Our results show 
that if the initial orbit of the outer binary A:(B+C) is circular (ei = 0), then the B:C pair 
is Hill stable for a wide range of e 2 (its initial eccentricity) e 2 < 0.37 and its semi-major 
axis a 2 < 26 AU. If ei = e 2 = 0, stability can be assures for < a 2 < 27 AU. Alterna- 
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tively, assuming cautiously that in the current configuration both orbits are at apocenters 
(maximum extension) we find that eccentricities e± < 0.19 and €2 < 0.20 are required for 
stability. In summary, the 3-body, {—EL 2 ) criterion guarantees that a nearly coplanar hi- 
erarchical system (with geometry consistent with VWCha), having small eccentricities, will 
never restructure by ejecting any star to infinity or form a close pair of either A:B or A:C. 

The analytical criterion cannot guarantee that the B and C components will not phys- 
ically collide. It might also indicate instability where none actually occurs on interesting, 
short, time scale (which we fix at 1 Myr). To confirm its predictions, we therefore integrated 
the motion of VW Cha-like configurations with a 7th-8th order Runge-Kutta method, accu- 
rate to 10 significant digits in the total energy. We assumed a\ = 110 AU and 0,2 = 16 AU. 
Only the ratio a 1 /a 2 ~ 7 : 1, is decisive for stability. Mutual approaches of stars to less than 
10 7 km (~ 10 PMS stellar radii) were treated as physical mergers following a tidal capture. 
We found that the analytical and numerical criteria approximately coincide. For instance, 
we confirmed numerically our analytical predictions discussed above. 

Fig. 2 presents numerical results for the planar and inclined system configurations. 
Curves are labeled with the starting initial inclination of the inner and outer orbits. Com- 
binations of the eccentricities located below the plotted curves represent stable systems. 
All destruction modes including physical collisions of stars were observed above the curves. 
Strictly speaking, the lines shown in the figure result from averaging of the results over such 
unessential initial parameters as the longitudes of pericenters, ascending nodes and true 
anomalies 

From Fig. 2 we conclude that a wide variety of eccentricities result in stable behavior. 
While ei < 0.45 is generally required, larger €2 values (up to 0.8) are consistent with stability, 
independent of the mutual orbital inclination (with exception of unstable, almost perpen- 
dicular orbits not presented here). We conclude that the hierachical projected structure of 
VW Cha likely reflects a stable physical configuration for the estimated age of the system. 

5. Conclusions 

The young double star VWCha (0"7 separation) has been successfully observed with 
adaptive optics techniques in the infrared J and K bands. These AO observations led to the 
followig main findings: 

• The secondary is a very close binary itself (O'.'l separation). 

• This physical triple system is shown to be stable over a period exceeding its estimated 
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age of 0.4 Myr. 

• Follow-up astrometric observations may, within a few years, constrain the dynamical 
mass of the B:C subsystem, thus testing evolutionary models of PMS stars. 

• The primary (A) has a significant IR-excess, indicating the presence of a circumprimary 
disk. This disk is probably small, with a radius not exceding ~30AU. 

Our VW Cha result shows that individual components of known binaries may be close 
binaries themselves, which needs to be considered while deriving ages and other parameters 
of such systems. 

We had many helpful discussions on VWCha with Goran Olofsson and Gosta Gahm. 
We are also grateful to Keith Hege for extensive support with the use of IDAC and to 
Wolfgang Brandner for supplying photometric data. 
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Fig. I. — An unprocessed frame in J (Left), and deconvolved images in J (Middle) and 
K (Right) of the VWCha system. The secondary itself is split into a 0.10" binary in the 
deconvolved images. The vertical bar in each image corresponds to one arcsecond. North is 
up, east is to the left. 
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Fig. 2. — Stability diagram of a model triple system representing VWCha. Eccentricities 
of a large outer and a small inner binary subsystem are denoted e±, e 2 , and their relative 
inclination is marked next to the curves. The system is stable at least 1 Myr for parameter 
sets below the plotted curves. 
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Table 1. IR Photometry and Astrometry of the VW Cha system 2000 03 11-13 



Stellar 


[J] 1.25 fim a 


[K] 2.18 fim h 


J-K 


Separation 


Position Angle c 


M bo i 


L 


M 


Age 


component 


(mag) 


(mag) 


(mag) 


(") 


(°) 


(mag) 


(L©) 


(M©) 


(Myr) 


A 


9.21 


7.03 


2.18 ±0.16 






3.56 


2.93 


1.00 


0.4 


B 


10.55 


8.92 


1.63 ±0.16 


0.661 ±0.006 


176.6 ±0.5 


4.97 


0.79 


0.40 


0.4 


C 


10.87 


9.32 


1.55 ±0.16 


0.100 ±0.006 


233.3 ±3.3 


5.29 


0.59 


0.35 


0.4 


B+C 


9.94 


8.35 


1.59 ±0.20 


0.683 ± 0.008 d 


179.4 ±0.7 d 


4.37 


1.38 


0.75 




A+B+C 


8.76 


6.75 


2.01 ±0.24 






3.14 


4.30 


1.75 





a Absolute error crj = 0.06 mag, mean relative error crj = 0.026 mag. 
b Absolute error ctk = 0.08 mag, mean relative error ctk = 0.12 mag. 

Reparation and position angle for components B and B+C are relative component A, and component C is relative component B. 
Position angle is measured from north to east. 

d For B+C the astrometry is made with respect to the photocenter. The larger error is due to the shift of the photocentcr with 
wavelength. 



